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We perform large shell model calculations for Calcium isotopes in the full fp shell 
by using the realistic Kuo-Brown interaction. The Calcium isotopes are especially 
interesting because the nearest-neighbour spacing distribution P{s) of low-lying 
energy levels shows significant deviations from the predictions of the Gaussian Or- 
thogonal Ensemble of random-matrix theory. This contrasts with other neighbour- 
ing nuclei which show fully chaotic spectral distributions. We study the chaotic 
behaviour as a function of the excitation energy. In addition, a clear signature 
of chaos suppression is obtained when the single-particle spacings are increased. 
In our opinion the relatively weak strength of the neutron-neutron interaction is 
unable to destroy the regular single— particle mean-field motion completely. In 
the neighbouring nuclei with both protons and neutrons in valence orbits, on the 
other hand, the stronger proton-neutron interaction would appear to be sufficient 
to destroy the regular mean— field motion. 

1 Introduction 

The nearest neighbour spacing distribution of energy levels provides a good 
signature of the chaoticity of a quantum system. In fact, the fluctuation prop- 
erties of quantum systems with underlying classical chaotic behaviour and 
time-reversal symmetry agree with the predictions of the Gaussian Orthogo- 
nal Ensemble (GOE) of random matrix theory, whereas quantum analogues of 
classically integrable systems display the features of the Poisson statistical. 

The nuclear shell model, with a realistic interaction and large configuration 
space, is one of the best theoretical approaches to the study of nuclear spectra. 
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The model provides large sets of exact energy levels and wave functions in 
truncated space, and their statistical analysis can give information on the 
features and borderlines of the transition from regular to chaotic dynamics 
in nuclei. The statistical analysis of shell-model energy spectra and wave 
functions has mainly concentrated on the sd shell region and chaotic behaviour 
has been found for these nuclei both near the yrast line and at higher energiefl. 

In this paper we extend xiur recent statistical analysis of the shell-model 
energy levels in the fp shelH, in order to study whether an order to chaos 
transition can be observed, and how it depends on variables such as excitation 
energy, angular momentum and single-particle energy spacings. 

2 Calculations in the fp shell 

The nuclear shell-model Hamiltonian, in second-quantization notation, can be 
written as 

iJ = ^ eafl+aa + < af3\V\6j > a^Opa^ag , (1) 

where the labels denote the accessible single-particle states, is the corre- 
sponding single-particle energy, and < Q;/3|y|(57 > is the two-body matrix 
element of the nuclear residual interaction. 

We follow the standard approach to obtain the eigenvalues of the shell- 
model Hamiltonian. Exact calculations for several nuclei are performed in the 
{f^/2^ P3/2: fb/2^ P1/2) Configuration space, assuming a ^°Ca inert core. The 
construction and diagonalization of large shell-model matrices are performed 
using a modified version of the computer code ANTOINEeI. For a fixed number 
of valence protons and neutrons we calculate the energy spectrum for total 
angular momentum J and total isospin T . The interaction is s, minimally 
modified Kuo-Brown realistic force with monopolc improvements. Coulomb 
effects are not included. 

To analyze the energy level fluctuations, it is necessary to consider only 
levels which have the same symmetries. In our case this means the same num- 
ber of nucleons and the same total angular momentum, parity and isospin. 
The level spectrum is then mapped onto unfolded levels with quasi-uniform 
level density. The suitable unfolding procedure depends on the region of the 
spectrum to be analyzed. For the low-lying levels, we use an unfolding proce- 
dure based on the constant temperature formula, where the mean level density 
can be assumed to be of the form 

p(i?) = iexp[(i?-i?o)/r], (2) 
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Table 1 : Brody parameter lo for the nearest neighbour level spacing distribution for < J < 
9, T = Tz states up to 4, 5 and 6 MeV above the yrast line in the analysed nuclei. 



Energy 


4bv+4t>Ti+*Sc 


4''Ca+^»Ca+='"Ca 


< 4 MeV 


0.92 


0.56 


< 5 MeV 


0.93 


0.60 


< 6 MeV 


0.95 


0.61 



with T and Eq fitting parameters. We liave compared this unfolding method 
with the standard local unfolding methocEI. Provided that the number of energy 
levels is not too small, the two procedures give similar results in the low-energy 
region, but the constant temperature level density is smoother and is preferable 
in the ground-state region. When the analysis includes many levels or the 
full spectrum, we use the local unfolding because, as is well known, due to 
the finite size of the shell-model basis, the eigenvalues are generally Gaussian 
distributed. 



3 Analysis of spectral statistics 

The P{s) distribution of the nearest-neighbour spacings Si = N [Ei+i) — N [Ei] 
of the unfolded levels is the best spectral statistics to study the fiuctuations of 
the short range correlations. To quantify the chaoticity of P{s) in terms of a 
parameter, we compare it to the Brody distribution, 

P(s,w) = a(w + l)s'^exp(-as"+^) , (3) 

with 

« = (r[^])-+^ . (4) 

This distribution interpolates between the Poisson distribution {lo = 0) of 
regular systems and the Wigner distribution {lo = 1) of chaotic ones (GOE). 
The parameter ld can be used as a simple quantitative measure of the degree 
of chaoticitja. 

We analyze the energy spectra looking at the lowest energy region, up 
to a few MeV above the yrast line. The main problem in this region is that 
the number of energy levels of the same symmetry is too small for a reliable 
statistical analysis. But having at our disposal the whole energy spectrum, we 
can use a sufficiently large set of levels to determine the secular behaviour of 
the level density and perform the unfolding procedure for each symmetry class 
of states. Moreover, in order to obtain more meaningful statistics, after the 



3 



1.1 




0.5 - 



0.4 I 1 1 1 < < 1 1 1 1 

4 6 8 10 12 14 16 18 20 22 

Energy (MeV) 

Figure 1: Energy dependence of the Brody parameter oj in **Ca. All levels up to an energy 
E above the yrast line are included for each J value. 

unfolding we can combine the level spacings of different J in a nucleus, or even 
of different nuclei, to calculate the P{s) distribution. 

A good estimate of the Brody parameter is obtained combining spac- 
ings of different nuclei. Table 1 shows the results for ^^V+^^Ti+'*^Sc and 
"^^Ca+^^Ca+^oCa, up to 6 MeV above the yrast lines. The number of level 
spacings is now sufficiently large to yield meaningful statistics and we see that 
Ca isotopes are not very chaotic at low energy, in contrast to other nuclei in 
the same region. 

The combined data of ^^Ca+^*Ca+^°Ca show that the chaoticity increases 
slightly with excitation energy up to 6 MeV. It is generally believed that states 
at higher energy, in the high density region, should be much more chaotic 
in nature. Thus, calcium isotopes offer the possibility to study the energy 
dependence of the chaoticity in the framework of the shell model. 

We have calculated the P{s) distribution and the Brody parameter up to a 
fixed value of the excitation energy above the yrast lines from 6 to 20 MeV. Wo 
use the local unfolding method and include all levels with J = 0-11. We find 
that the chaoticity increases rather smoothly with energy. This is illustrated 
in Figure 1 for ''^Ca. It is necessary to include levels up to about 14 MeV to 
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Figure 2: P{s) distribution for **Ca. All levels for each J value are included. The single- 
particle energies axe: degenerate, normal and double spaced. 



get Lu > 0.9. 

Table 2 shows the Brody parameter u) for the whole spectrum of the anal- 
ysed Ca isotopes, which range from ^^Ca to ^"Ca. We see that the lightest Ca 
isotopes are not fully chaotic even when the whole energy spectrum is taken 
into account. 

For the heavier calcium isotopes the number of states is very large, e.g. 
17,276 for ^°Ca. Therefore, it is possible to analyze separately the spectra 
for different J values with good statistics. We do not find any significant 
dependence on angular momentum, except for J = where we find w ~ 0.8 
both for "'^Ca and ^°Ca. This is probably related to the fact that the pairing 
interaction, which preserves the seniority quantum number, is more effective 
for the J = states. 

It is very interesting to study the effect of the one body Hamiltonian on 
the P{s) distribution. The single-particle motion in the spherical mean field is 
regular, while the nuclear two-body residual interaction is strongly non-linear. 
Figure 2 shows how the fluctuation properties of nuclear energy levels in ^^Ca 
change when the single-particle spacings are changed. We consider three cases: 
i) the /7/2, P3/2) /5/2 and pi/2 single-particle spacings are all degenerate; ii) 
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Table 2; Brody parameter w for the nearest neighbour level spacing distribution for Ca 
isotopes. All levels for each J value are included. 



44Ca 


"t-Ca 


*Ca 


*^Ca 


4«Ca 


«Ca 


''^'Ca 


0.69 


0.75 


0.99 


0.98 


0.95 


1.00 


0.87 



the single-particle spacings are the experimental ones; iii) the single-particle 
spacings are multiplied by a factor of two. The realistic residual interaction 
is the same in all cases. The P(s) distribution includes all the nuclear levels, 
separately unfolded for each J value. For degenerate single-particle levels 
Figure 2 clearly shows chaotic dynamics in ^^Ca, but when the single-particle 
spacings are increased, there is a transition towards regularity. The same kind 
of chaos to order transition is observed for the other Calcium isotopes. 

4 Conclusions 

The Ca isotopes are especially interesting because the nearest-neighbour spac- 
ing distribution P{s) of low-lying energy levels shows significant deviations 
from GOE predictions, in contrast to other neighbouring nuclei which show 
fully chaotic spectral distributions. 

The analysis of level spacings up to a given excitation energy shows that 
the chaoticity of Ca isotopes increases smoothly with the excitation energy end 
point. However, even when the whole energy spectrum is included, the lighter 
isotopes ^^Ca and "^^Ca are not fully chaotic, and the order to chaos transition 
is progressive as the number of active particles increases. 

In the heavier Ca isotopes studied, the number of energy levels is very large 
and it is possible to analyse with good statistics the energy levels of a single 
J value. For example, in ^"Ca the dimension of many Hamiltonian matrices 
exceeds one or two thousand for fixed J. We do not find significant differences 
in the P{s) distribution for different J values, except perhaps for J = 0, for 
which we obtain a Brody parameter somewhat smaller than for other J values. 
This is probably related to the pairing component of the interaction and the 
aproximate conservation of the seniority quantum number. 

Finally, a clear chaos to order transition is observed as the single particle 
spacings are increased. In the lighter Ca isotopes we obtain chaotic behaviour 
for degenerate single-particle states and quasi-regular motion when the single 
particle spacings are twice the experimental values. Thus, it seems that the 
main reason for the substantial deviations from chaoticity obtained in Ca iso- 
topes is that the relatively weak strength of the neutron-neutron interaction is 
not able to destroy the regular single-particle mean-field motion completely. 
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But in the neighbouring nuclei with both protons and neutrons in valence or- 
bits, the stronger proton- neutron interaction seems to be sufRcent to destroy 
the regular mean-field motion. 
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